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Hammersley’s  Law  for  the  van  der  Corput  Sequence: 

An  Instance  of  Probability  Theory  for  Pseudo -Random  Numbers 

by 

A.  del  Junco  and  J«  i'iichael  Steele 


1.  Introduction, 

Pseudo-random  numbers  form  the  backbone  of  computer  simulation  and 
Monte  Carlo  analysis,  yet  there  are  essentially  no  known  theorems  which 
make  explicit  the  sense  in  vftiich  pseudo-random  ntimbers  are  replacements 
for  random  numbers.  There  is  probably  nq  general  result  which  can  be 
proved,  but  there  is  still  a program  by  -which  a deeper  xinderstanding  of 
the  relationship  between  pseudo-random  numbers  and  random  nuiribers  can 
evolve. 

For  every  pseudo-random  sequence  and  every  probabalistic  theorem^ 
there  is  surely  some  analogue  of  the  probabalistic  -theorem  for  -that 
pseudo -random  sequence.  The  qualities  of  -that  analogue  sho-uld  then 
accurately  reflect  the  random  qualities  of  the  pseudo-random  sequence. 
Ry  a systematic  analysis  of  -these  pairs  of  theorems  and  sequences,  a 
body  of  results  can  be  obtained  which  is  capable  of  resol-ving,  or  at 
least  eroding,  many  of  the  philosophical  and  practical  questions 
concerning  pseudo-random  sequences. 

This  is  no  over-night  task  and  there  is  no  ob-vious  place  to  begin 
other  -than  -wi-th  the  theorems  and  sequences  -vdiich  interest  one  most.  ¥e 
have  begun  with  a -theorem  due  to  Hammersley  and  a sequence  due  to 
van  der  Corput. 
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be  indepen- 


To  introduce  Hammersley' s theorem,  let  X^,  i = 1,2,..., 
dent  random  variables  vith  xiniform  distribution  on  [0,1],  "By 

we  denote  the  cardinality  of  the  largest  monotone 
increasing  subsequence  of  the  values  Xj^(a3),X2(c3),,.,,X^(aj).  Hammersley 
[2]  proved  that 

(l.l)  lim  = c 

^ n 

n oo 

where  c is  a constant  and.  the  convergence  is  in  probability.  Vie  axe 
particularly  interested  in  (l.l)  because  of  the  considerable  effort  xdiich 
has  been  focused  on  the  determination  of  c.  Even  before  convergence  had 
been  proved  in  (l.l),  Baer  and  Brock  [l]  had  conjectured  that  c = 2 on 
the  basis -of  extensive  computation.  Hammersley  [2]  gave  bounds  on  c 
which  were  improved  by  Kingman  [4],  but  the  deepest  results  obtained  so 
far  are  due  to  Logan  and  Shepp  [6]  who  proved  c > 2,  The  problem  of 
proving  c < 2 remains  open. 

In  the  present  context  Hammersley’ s theorem  has  a natural  appeal 
as  a substantive  probabalistic  result  in  which  there  is  much  current 
interest.  The  tempting  prospects  of  obtaining  c in  the  analogue  of 
(l.l)  for  some  pseudo-random  sequence  provide  even  more  reason  to  start 
with  Hammersley’ s theorem. 

One  natural  choice  for  the  sequence  of  pseudo-random  numbers  would 
be  those  gen.erated  by  the  linear  congruence  method.  These  are  the  most 
widely  used  pseudo-random  numbers,  yet  they  are  very  difficult  to  analyse. 
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For  this  reason,  -we  have  begun  vri.th  the  van  der  Corput  sequence  vftiich  has 
•wide-spread  use  in  numerical  integration  [3]  as  -well  as  considerable 
historical  interest.  (See,  for  example,  Knu-th's  interesting  discussion 
"’ttiat  is  a random  sequence"  [^5,  pp.  127-157]  )» 

To  define  -fche  van  der  Corput  sequence  -we  first  write,  for  n > 0, 
n = i=0  "Where  a^  = 0 or  1.  The  n+lst  element  of  "the  sequence 

is  ^>2(11)  = X i=o  • One  can  see  that  ==  fPgCl)  = ^ / 

1 3 

= 5"  ^ The  nature  of  cp2(j“)  is  more  easily  seen 

in  binary  notation  -vAiere  i®  "the  refection  of  n in  the  decimal 

point."  As  the  subscript  suggests  one  can  define  a similar  sequence  by 
= X i=o  where  n = ^ a^p^,  0 < a^j^  < p and  P > 2 is 

an  integer. 

We  no"W  have  the  theorem: 

If  f (n)  is  the  cardinal! -ty  of  the  largest  monotone  increasing 

subsequence  of  {9  (O),  cp  (l),...,(p  (n-l))  then 

P P P 

(1.2)  lim  n~^^^  ^ (n)  t lim  3a  -?(n)  = ^ for  p = 2 

n ->■  oo  n *+■  00 


and 

(1.3)  11m  n~^^^  ^ (3a)  =2  -/l-p  ^ , lim  n I (n)  = p^^^  for  p > 2 . 

n -*  CO  n “►  CO 

This  result  is  as  precise  an  analogue  to  (l.l)  as  one  could  realistically 
expect.  Also,  the  Logan-Shepp  loiver  bound  on  c masses  it  particiilarly 
notetrorthy  that 
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(lA) 


i (n)  = 2 . 


lim  lim 

p CO  n 00 

For  digestibility,  the  proof  of  these  results  has  been  divided  into  fotir 
parts.  First  we  obtain  a geometrical  characterization  of  the  monotone 
subsequences  of  van  der  Corput’ s sequence.  Upper  and  lower  bounds  are 
then  obtained  for  i (n).  Finally  the  required  limits  are  identified. 

2.  The  Geometrical  Representation. 

Our  first  goal  is  to  obtain  a representation  from  \jhxdh  ^tailed 
information  about  C<Pp(0 ),9^(l), . . .,cp^(n-l))  can  be  deduced.  We 
define  cr^  to  be  the  unique  permutation  of  [0,1, . . . ,n-l}  which  puts 
9 (i),  0 < i < ii*’!  in.  increasing  order,  i.e. 

'P  (o'  (o))  < *P  (o'  (l))  < ...  < 9 (cr  (n-l))  . 

This  permutation  mil  be  written  as  a sequence 

0-^  = ■(o'jj(0),o-^(l),...,o-j^(n-l)) 

and  our  reason  for  introducing  cr^  . is  the  elementary  fact  that  the  length 

of  the  longest  monotone  increasing  subsequence  of  cr^  is  equal  to  f (n). 

We  also  have  the  following  law  of  formation  of , cr  ^ vhich  will  be 

mp 


crucial. 


Leama  2.1. 


One  obtains  cr  ^ from  o'  replacing  the  entry  cr  (i) 

mp  p 

of  cr  by  the  sequence  cr  (i)+p  ct  . 
n n m 

P P 

Remark.  Here  and  in  the  following  if  x is  a sequence  then  a-^bx  is  a 
sequence  with  the  same  domain  as  x defined  by  (a+bx)(i)  = a+bx(i). 

Proof  of  lemma.  Notice  that  there  are  two  things  we  must  show: 


(2-1)  9p(cr  ^(i)  + P^o-^Ci))  < 9p(cr  j^Ci)  + p’^cr^(d-HL))  , 


and 


(2-2)  cp  (a  ^(i)  + p”a-^(d))  < 9 (o’  n(i+l)  + , 


for  any  j,  j*. 

We  first  suppose  j < p”  and  k is  any  integer.  Setting  d = T a-P^^ 

^ 1=0  1 

^ _ y m a^p^  ^ we  have  j + p\  = ^ a^p^,  and  consequently 

~ ^ i=n 


(213) 


9 

P i=0 


-(i+1) 


4.  V “(i+l) 

= ^a.p  +P  2l®-P  ' 

i=0  ^ i=n  ^ 


= 'Pp(o)  + p”%p(k)  for  j < p^ 


Now  by  (2.3)  we  see  (2.1)  is  equivalent  to 


V~\^(%(3))  < P~“9p(o-^(D+l))  , 


-n. 
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■'.chich  is  just  the  dsfiniiion  of  cr^.  Siailarly,  (2.2)  is  equivalent  to 


(2A)  cpp(cr  ^(j))  + p'^^cPpCo-^Cj))  < <Pp(cr  (o-^^Co'  )) 

P JP 


To  Chech  {2,h)  we  note  that,  as.  i runs  throu^  {0,1, p^-l},  q>p(i) 


runs  throu^  Co^P  ^,2p"”', (p^-l)p  “}  so  that,  in  fact, 
cp  (cr  ^(i+l))  = p ^ + 9^(o'  ^(i)).  The  proof  of  the  lemma  is  thus  complete. 

In  many  ways  it  is  easier  to  work  with  the  permutation  matrix 
associated  with  cr^.  We  recall  that  is  an  nxn  matrix  defined 

by 


n , N.  -n-, 


if  i 


0 otherwise  . 


To  rephrase  Lemma  2.1  in  terms  of  matrices  we  define  m different 

p’^  X mp^  matrices  ^ 'by 

P^  P P 


and 


P 


( A j^(i,jm  if  m(  j 

’ 

1 otherwise 


A^^(i,j)  = A*^  (i,  (j-i)mod  m)  . 
P P 
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Intuitively,  we  get  by  inserting  m-l  columns  of  O's  after  each 

P -0 

colnmn  of  A ^ and  then  shift  A ^ by  ^ places  to  the  right  to  get 


T> 


P 


A' 


n 


Finally  we  can  jhrase  Lemma  2.1  by  the  following  formula  for  A 


as  a block  matrix: 


mp 


(2.5) 


A 


n 


mp 


P 

o-“^(2) 

A 

n 

P 


P 


This  matzux  representation  of  the  permutation  permits  a geometrical 

view  of  the  decreasing  subsequences.  Formally  an  increasing  subsequence 

of  (T  is  a sequence  of  integers  P = <i„}  such  that 

n X d.  ^ 

) < o'(ji+p)*  Such  a subsequence  can  be  identified  in  A^  as  a path 
of  I's  which  goes  down  and  to  the  ri^t.  The  set  of  integers 

= (k:  ^ < 3^)  will  be  called  the  domain  of  the  path  P 

and  will  be  denoted  by  dom(p).  The  integer  interval  [cr  (3-,  )^cr  (3,)]  = 
{k:  cr  (3n  )<k<cr(3.)}  will  be  called  the  range  of  the  path,  and 
naturally  i (P)  will  be  called  the  length  of  P. 
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3.  The  lover  Bouad  lemma. 


The  main  result  of  this  section  is  the  foUoi-iing, 
lemma  3.I.  For  ^ < m < ^ ve  have 

(5.1)  ^ (mp’^)  > + m(p-l)  - (3p-l)  , 

, ^ n-3  ^ ^ n-2 

and  for  p < m < p ve  have 


(3.2) 


;8-(mp'^)  > p^~^(p-l)  + mp  - (p^+2p-l) 


Proof.  First  ve  vill  consider  m such  that  p^~^  < ^ 

XI  IL. 

construct  an  explicit  path  throu^  A By  setting  lo-  = (p  “ -m)/(p-l) 

2 

and  ” ^“®i  “have  m^p  + m^p  = p^.  For  0 < i < [m^]  we  can  define 

2 2 

a path  P.  throu^  A of  length  2p-l  and  domain  [p  i,p  (i+l)-!]  by 

P 

P.  = (p^i,p^i-Jp,p^i+2p,...,p^(i+l)-p,p^(i+l)-p+l,...,p^(i+l)-l)  . 

X 


For  0 < i < [n^]  we  then  define  a path  Pr  thro^3^  A ^ with 

2 2 ^ 
length  p and  domain  [[m^]p  +ip, [m^]p  +(i+l)p-l]  by 


Ifow  we  consider  the  paths  P.  in  A ^ which  correspond  to  P. . 

Fnese  paths  are  defined  by  P^  = mP^+cr^^(i).  Finally  we  define  our 
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desired  patli  P throng  A as  the  concatenation  of  the  P , 

n i'' 

mp 

P = P_Pt  ...  Pr  T.r  It* 

0 1 

To  make  certain  that  the  path  P is  properly  defined  we  have  to 
check  that  for  each  k e dom  P.  and  £ e dom  P.  we  have  k > Ihis 
condition  will  be  abreviated  by  dom  P.  > dom  P. , and  it  follows 

1+1  i'^ 

immediately  from  the  corresponding  fact  that  dom  ^ ^ 

construction  we  finally  have 

^(mp’^)  > je(P)  = [m^](2p-l)  + [B^]p 

> n^(2p-l)  + m^p  - (3p-l) 

= p”~^  + m(p-l)  - (3p-l) 


which  completes  the  proof  of  (3.1). 

We  now  consider  m^  p^  ^ < m < p^  ^ and  begin  by  setting 

= (p^”^-ni)/ (p-l)  and  5^  = m-m^.  This  time  we  have  m^p^-hr^p^  = p^. 

2 

As  in  our  previous  construction  there  is  a path  of  length  p through 

2 

A so  for  0 < i < [m^]  there  is  a path  P^  of  length  p throu^ 

!P  Zt  "T 

A with  domain  [ip^ , (i+l)p^-l].  For  0 < i < [mu]  there  is  a path 

pU  d. 

Pr  T of  length  2p-l  throu^  A mth  domain 
I m,  j +x  n 

X 2 3 2 ^ — —1 

[[m^ljK+ip  ^ [in^]p'^+(i+l)p  -1].  Next  we  set  P^  = mP^  + cr~  (i)  and 

O’  ^(i) 

— m ^ 

obtain  a path  throu^  A by  letting 
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P — 1 if  In* 

0 1 [ns^J+lBigj-l 

Finally  we  calculate  as  before: 

Z (mp^)  > i (P)  = (2p-l) 

> iii^p^*Hi^(2p-l)-(p^+2p-l)  = p^"^Cp-l)+n^-(p^+2p-l) 
which  conroletes  the  proof  of  (5-2) . 

Remarh.  The  preceeding  proof  probably  appears  more  iamalved  than  an 

examination  of  an  example  would  indicate.  Ey  calctilating  A for 

mp 

p = 3 and  several  small  values  of  m and  n one  can  easily  find  the 

path  P and  see  how  it  evolves  as  m and  n increases.  Similarly  in 

the  arguments  idiich  follow  one  should  keep  the  example  p = 5 clearly 

in  mind,  perhaps  by  keeping  a small  table  of  the  A 

mp 

h.  The  Upper  Bound  Lemma. 

The  recursive  nature  of  A provides  the  key  to  the  following, 

mp“ 

Lemma  4,1.  For  any  positive  integer  m we  have 
(4.1)  Z (mp^"^)  < p^~^  + m(p-l) 


and 


(4.2) 


. , nn  ^ n-1  n-2  , 

Z (mp  ) < p - p + mp  . 
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we  can  decompose  P as  a 


Proof.  For  any  path  P throii^  A 
concatenation 

where  P^  is  a path  throi:^ 
range  of  P^  is  a subset  of 

We  now  let  P.  be  the  corresponding  path  in  A of  P. . That  is, 

pH  X 

we  let  P.  = m ^ (P.  -cr  ^ (i ) ) - We  do  not  have  dom  P.  < dom  P.  ^ but  we 

obviously  do  have  dom  P^  < dom  so  the  two  paths  have  at  most  one 

common  point.  If  dom  P.  H dom  P.  ^ we  shall  say  P.  and  P.  are 

1 1+1  ' ^1  1+1 

El  1 

linked.  Since  f (P)  < -2  (P^^)  the  lemma  would  follow  from 

^ i^  f (Fj_':i)  < min(p°“^-ha(p-l),p”“^-p”~^-ffflp).  This  fact  will  be  proved 

using  only  the  inequality  dom  P^  < dom 

First  we  decompose  A into  blocks  of  hei^t  and  width  p. 

pn 

We  note  that  each  block  bas  exactly  one  1,  and  the  pattern  of  l*s  in  each 
column  of  blocks  is  the  same  as  the  pattern  in  Ap. 

The  block  width  w(P)  of  a path  P is  defined  as  the  number  of 
distinct  columns  of  blocks  -^diich  intersect  P.  The  block  hei^t  h(P) 
is  defined  correspondingly.  Further,  we  let  p equal  the  number  of 
integers  i such  that  P.  is  linked  to  P.  , , . 

If  P.  and  P^  are  linked  we  define  their  linkage  as  the 

concatenation-  P. Pf,,  where  Pf,,  denotes  P.  -v/ith  its  first  entry 
removed.  By  forming  the  successive  linkage  of  the  P^  vre  obatin  m-p. 

A A A 

paths  P,  ,P^,...,P  which  are  unlinked  and  for  which 
1 ' m-u 


mp 


P = P P 
^0  1 


,.  P 


m-1 


7 

A 

n 
P 

Cip° ,(i+l)p^-l], 


Here  one  should  note  that  the 
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(^.3) 


^ m-l 
i=l  i=0 

Next  we  fom  all  possible  concatenations  of  successive  R to  get  s paths 
. . jPg  ■viiida  allow  no  further  concatenation.  One  then  has  s < m-p 

and 

s m-p  . 

5;  i(Pi)  = £ i(P.)  . 

1=1  i=l 

MDreover,  one  notes  that  no  two  P^  share  a conmon  column  of  bloclcs  for 
the  simple  reason  that  dom  P^  < dom  ^ ariy  block  is  to 

the  ri^t  of  the  1 in  any  block  above  it. 

The  crucial  observation  is  that  for  each  path  P^  we  have 

j^(Pp  <h(Pp  +w(Pp-l  <p  + w(Pp-l 

viiich  by  summing  over  i yields 

(^.5)  Z^(Pi)  < E w(P')  + s(p-l)  . 

i=l  i=l 

Since  no  tiv’o  P^  share  a column  of  blocks  we  have 

(i^.6)  ^wCpp^p"^"^. 

i=l 

By  the  inequality  (^.5)  obtain 

(.h,7)  Z ^(P' ) < p""^  + s(p-l) 

i=l 
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■which  h-y  setting  s = n-ii-X  -with  X > 0 -using  and  (h,k-)  hecomes 


(h,8)  X (m-X)(p-l)  - h(p-2)  . 

i=0 

Ttiis  completes  the  proof  of  (if-.l)  since  X > 0 and  h > 0.  On  the  other 
hand  (^.6)  also  yields  s < p^  so  -JJ-  < p^^~^-(m-X)  -which  by  (^.S) 

implies 


i ^(P  ) < p““^  + (m-X)(p-l)  - (p’^'^-(m-X))(p-2) 
i=0 


. n n-1 , 

< p -p  -ha 


This  iurolies  £(m>^)  < p^-p^'"^-ita  and  by  replacing  n by  n-1  and  m 
by  mp  ve  obtain  the  proof  of  (^.2).  This  completes  our  proof  of  the 
upper  bound  lemma. 

Before  appiying  the  preceeding  lemmas,  -we  should  note  that  in  certain 
cases  the  inequalities  can  be  made  to  pro-vide  equality.  By  inspecting 
the  proof  of  (5.1)  in  the  cases  of  m = p^  ^ and  m = p^~^  one  sees 
that  it  is  not  necessary  to  ^estimate  greatest  integers  so  that  the  3P“1 
disappears  from  the  fi^t  of  (3.1).  Combining  this  -svith  (^i-.l)  i<re  have 
■the  identities 


(^•9)  i(p^“"^)  - P^  and  = p^~^(2p-l)  . 

By  similar  analysis  for  p = 2 ve  also  have  the  identity 
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(U.IO) 


£(e2^) 


for 


,n-2 


< m < 2 


n-1 


2”“^+2m 


for 


.n-3 


< m<  2 


n-2 


5'.  Identification  of  the  T.imi-hH. 

To  conrolete  the  proof  of  our  theorem,  we  first  identify 

where  N is  restricted  to  S = {115)”:  ^ £ p^~^}. 

Lemma.  5*1*  Setting  K = max(2(l-p~-^)^'^^,  2-p~^)  and 

k = min(2(l-p  2-p  we  have 


lim  (F) 

n -+  CO 


(5.1) 


lim  jg  (F)  = K 

FeS 


and 

(5.2)  ^ f"^/^  f (F)  = k 

FeS 

Proof.  As  in  Section  3 w-e  deal  with  p”  ^ < m < p^”^  (Case  l)  and 
p < m < p^~  (Case  2)  separately.  This  time  we  begin  with  case  2. 
% (3.2)  and  (^.2)  we  have 

gjj(ia)-(p^+2p+l)(mp")''^^^  < jg(mp^)(mp^)"^/^  < 

n-?  n 

where  = (p  (p-l)+xp)  (xp  ) 


lii- 


. Consequently  we  have. 


(5.3) 


lim  2,  (mp^)(mp^) 


-1/2 


= 3^ 


where  the  lim  inf  is  tahen  over  the  set  S = < 


m 


< n3 , 


One  haus  the  same  equation  as  (5 .5)  for  the  3im  sup. 

Next  observe  has  its  mirdmum  at  = p^~^(p-l)  anr?  is 

decreasing  for  p^”^  < x < a and  increasing  for  a <a:<  One  easily 

— — n n — — 

checks  that  g„(p““^)  = p^^^,  g„(p“'^)  = 2-p"^,  and  g„(a„)  = 2Cl~p"^)^/^ 

il  Xi. 


■'n' 


so  we  have 


(5.M 


lim  g (m)  = k and  lim  g (m)  = K 

n „ n ^ 

mt>  eS, 


mp 


To  prove  the  comparable  identities  for  case  1 one  sets 

-n  “I/2 

fj^(x)  = (p  ■^+x(p-l))(xp  ) and  obtains  from  (5.1)  and  (^.1)  that 


lim  i(mp^)(mp’^)  = lim  f (m) 


mp^eS, 


mp  eSg 


n 


where  : P^  < m p^  n > l].  The  same  equality  holds  for  the 

limsupj  and  we  note  as  before  that  f^  has  a minimum  at  p ” (p-l)"  = b^^, 

fecreases  for  P°  ^ < x < b^  and  increases  for  b^  < x < p^"^.  Finally 
we  note  ?n(\)  = 2 (l-p“^)^'^^,  fj^(p^~^)  = and  f^(p“"^)  = 2-p"\  . 

Hence^  we  have 


(5.5) 


lam  I (m)  > k and 
mp  sS 


lim  f (m)  < K 
n - 

nip 
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One  has  inequality  in  (5*5)  since  is  not  an  integer.  Actually 

equality  can  be  proved,  but  is  not  required  for  the  rest  of  the  proof. 

^ (5-3 ),  (5.^)  and  (5.5)  we  have  completed  the  proof  of  (5.1)  and  (5-2) 
as  required  by  the  lemma. 

The  proof  of  our  theorem  can  now  be  completed  in  a routine  way.  For 
an  arbitrary  integer  N we  write 

rc  -V-  -4.-U  ’^“3  ^ ^ n-1  ^ ^ ^ n+2 

(5-6)  h = E5>  +r  with  p ^ < P 0 < r < p 

Letting  .2(i,o)  denote  the  cardinality  of  the  largest  increasing 
subsequence  of  (tp  (i);  9 (i+l), ...,<p  (j-1))  we  see 

Jr  Jr  Jr 

(5-7 ) (Ip^ , (i+l)p^ ) = -S  (p^) 

as  a consequence  of  (2.3).  There  is  also  the  obvious  fact  that  i(i,3) 
is  subadditive. 


Z (i,k)  < Z (i,d  ) + jg  (l^h)  for  i < j < k . 


n 

mp 


choosing  an  integer  t such  that  tp^  < mp^ 
+ p^'~  < (t-J-2)p^^^  one  has  by  (5-6)  that 


and 


l6 


s,  (K)  < Z 

< Z (mp  ')  + Z (fflp”^Jiip^+p”^^) 

< S,  (t+2}p^'^^  ) 

< Z (mp”)  + 2Z  (p”"^^)  . 

Hence  we  have, 

^5-8)  Z (mp”)  < i (H)  < i (mp”)  + 2Z  (p^"^^)  . 

Since  for  large  N one  has  N/mp”  is  near  1 and  f (p^'^^)/(mp”)^/^ 

IS  near  0,  the  inequality  (5.8)  together  with  Lemma  5 completes  the 
proof  of  the  theorem. 


17 


References 


[1]  Baer,  R.M.  and  Brock,  P.  (19^8),  "Natural  sorting  over  permutation 

spaces,"  Math.  Comp.  22,  385-^10* 

[2]  Hammersley,  J.Ii.  (197’2),  "A  few  seedlings  of  research,"  Proc.  Sixth 

Berk.  Symp.  Math.  Stat.  and  Proh.  ^ 3^5-59^. 

[3]  Hammer sley,  J.M.  and  Handscomh,  D.C.  Monte  Carlo  Methods,  Methuen, 

london, 

[4]  Kingman,  J.F.C.  (1973)^  "Subadditive  ergodic  theory,"  Ann,  Prob.  ^ 

883-899. 

[5]  Knuth,  D.E.  ihe  Art  of  Computer  Programming,  Vol.  2,  ( S emi numerical 

Algorithms),  Addison-Wesley,  Palo  ALto,  I989. 

[6]  Logan,  B.F.  and  Shepp,  L.A.  (1975)^  ”A  variational  problem  for 

random  Young  Tableaux."  Special  Invited  Paper  Annual 
Meeting  of  E-B.  To  appear  in  Adv.  Math. 


18 


UNCLASSIFIED 


security  CLASS! PICATION  of  THiS  PAGE  Data  Enttrad) 


READ  INSTRUCTIONS  I 

BEFORE  COMPLETING  FORM 

1 1.  REPORT  NUMBER 

25 

2.  GOVT  ACCESSION  NO, 

3.  RECIPIENT'S  CATALOG  NUH8ER 

} 4.  TITLE  (and  Subtitle) 

j Hammersley's  Law  for  the  van  der  Corput  Sequence: 
An  Instance  of  Probability  Theory  for  Pseudo- 
Random  Numbers 

5.  TYPE  OF  REPORT  a PERIOD  COVERED 

TECHNICAL  REPORT 

S.  PERFORMING  ORG.  REPORT  NUMBER 

7.  AUTHORfa; 

A.  del  Junco  and  Michael  Steele 

».  CONTRACT  OR  GRANT  NUMBERfaJ 

DAAG29-77-G-0031 

9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Department  of  Statistics 

1 Stanford  University 
j Stanford,  CA  9^305 

10.  PROGRAM  element.  PROJECT,  TASK 
AREA  a WORK  UNIT  NUMBERS 

P-114-14.35-M 

|l1.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

U,  S.  Army  Research  Office 

Post  Office  Box  12211 

Research  Triangle  Park,  NC  27709 

12.  REPORT  DATE 

October  9,  1978 

13.  NUMBER  OF  PAGES 

18 

U,  MONiTOHtHD  AGENCY  NAME  & ADDRE5S('i/  dlltBrent  from  Contrelting  OIIIcb) 

IS.  SECURITY  CLASS,  (ol  tMe  report) 

UNCLASSIFIED 

is.  DISTRiaUTlO.N  STATEMENT  (ot  Mm  Kapott) 


Approved  for  Public  Release^  Distribution  Unlimited. 


17.  DISTRIBUTION  STATEMENT  (ot  lha  abatract  antartd  In  Stock  30,  It  dlHatani  horn  Kaport) 


10.  SUPPLEMENTARY  NOTES 

The  findings  in  this  report  are  not  to  be  construed  as  an  official  Department 
of  the  Army  position,  unless  so  designated  by  other  authorized  docviments. 

This  report  partially  supported  under  Office  of  Naval  Research  Contract 
NOOOl^ -76 -C -0^4-73  (NR-042-267  ) and  issued  as  Technical  Report  No.  263.  

KEY  WORDS  (Contlnv0  on  sld*  It  Ti«cmB99fy  vtjd  identity  by  block  ttvtnbBr) 

van  der  Corput  sequence,  monotonic  subsequence. 


20.  ABSTRACT  (Continue  on  tevreo  elde  U necoBBory  end  identify  by  block  m2int>#f> 

Please  see  reverse  side. 


1473  edition  OF  1 NOV  9»  IS  OBSOLETE 

S/N  0102-014-6601  i 


nn 

UU  1 j an  73 


UNCLASSIFIED 


The  analogue  of  Hammersley' s theorem  on  the  length  of  the  longest 
monotonie  subsequence  of  i.i.d.  continously  distributed  random 
variables  is  obtained  for  the  pseudo-random  van  der  Corput  sequence. 
In  this  case  there  is  no  limit  but  the  precise  liminf  and  limsup 
are  determined.  The  constants  obtained  are  closely  related  to  those 
established  in  the  independent  case  by  Logan  and  Shepp,  and  Vershik 
and  Kerov. 


UNCLASSIFIED 


SECUKITV  CLA$tlFICATIOM  OP  THIS  PAOefiniMi  Data 


